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1  SUMMARY 


The  goal  of  this  research  project  is  to  develop  new  advanced  numerical  methods  and  to  per¬ 
form  DNS  studies  of  transient  hypersonic  reacting  flows  over  full  3-D  maneuvering  vehicles. 
The  DNS  tools  and  supporting  theoretical  approaches  are  used  to  gain  new  fundamental 
understanding  of  transition  phenomena  of  3-D  chemically-reacting  hypersonic  boundary  lay¬ 
ers.  In  the  three-year  period,  significant  progress  has  been  made  in  the  DNS  of  hypersonic 
l)oundary-layer  flow  instabilities  in  the  following  areas:  1)  receptivity  of  3-D  hypersonic  flows 
over  blunt  elliptical  cones,  2)  nonlinear  breakdown  of  3-D  Gortler  vortices,  and  3)  reacting 
hypersonic  flow  receptivities  and  stabilities.  We  have  also  finished  the  development  of  linear 
stability  codes  for  hypersonic  flow  and  have  started  the  work  on  developing  new  parallel 
implicit  algorithms  for  DNS  of  reacting  hypersonic  flows.  So  far,  we  have  developed  new 
iiK'thodologics  and  have  demonstrated  the  feasibility  of  our  approach  as  a  practical  DNS 
tool  to  explore  transient  hypersonic  reacting  flow  physics.  The  focus  of  future  work  is  on  the 
I).\S  studies  of  the  hypersonic  flow  physics  related  to  the  hypersonic  boundary  layer  transi- 
t  ion.  and  to  extend  these  DNS  and  theoretical  studies  to  3-D  transient  hypersonic  flows  with 
complex  real-gas  models  and  complex  geometries,  and  the  DNS  of  nonlinear  breakdown  and 
interactions  of  instability  modes  in  hypersonic  boundary  layers. 
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2  OBJECTIVES 


Objectives  and  Tasks 

Hypersonic  boundary  layer  laminar-turbulent  transition  and  unsteady  hypersonic  aero¬ 
dynamics  are  fundamental  problems  which  have  important  practical  applications  m  design 
and  control  of  hypersonic  lifting  vehicles  but  are  currently  not  well  understood.  In  recent 
years,  direct  numerical  simulation  (DNS)  has  become  a  powerful  tool  in  the  studies  of  the 
.stability  and  transition  of  boundary  layers.  However,  available  DNS  methods  cannot  be 
ajiplied  to  hypersonic  boundary  layers  over  realistic  blunt  bodies  because  the  effects  of  nose 
blunt  ness,  the  presence  of  shock  waves,  and  the  real-gas  effects  at  high  temperatures  are 
ncgh'cted  in  the  conventional  DNS  approach. 

Tli('  objective  of  this  research  project  is  to  develop  new  numerical  methods  and  perform 
D.XS  studies  of  laminar-turbulent  transition  of  hypersonic  boundary  layers  and  transient 
li>  personic  flows  over  complex  3-D  maneuvering  vehicles.  The  DNS  numerical  tools,  as  well  as 
ot  her  t  heoretical  approaches  such  as  the  linear  stability  analysis  and  the  parabolized  stability 
(■(|uations.  are  used  to  gain  a  new  fundamental  understanding  of  transition  phenomena  of 
3-1)  chemicallv-reacting  hypersonic  boundary  layers.  In  order  to  reach  the  research  goal,  the 
fulhnviug  tasks  were  undertaken  during  the  report  period: 

1.  To  develop  and  validate  new  high-order  massively-parallelized  numerical  methods 
and  parallel  computer  codes  for  DNS  of  general  3-D  hypersonic  reacting  flows  over 
realistic  maneuvering  hypersonic  vehicles.  The  new  methods  include  new  parallel 
methods  for  time-accurate  v'iscous  flow  simulations,  high-order  spatial  schemes, 
third-order  semi-implicit  Runge-Kutta  schemes  for  reacting  flow  equations,  and 
semi-implicit  spatial  schemes  for  viscous  flow  equations. 

2.  l  o  perform  DNS  studies  of  the  receptivity  and  stability  of  3-D  hypersonic  bound¬ 
ary  layers  over  elliptical  blunt  cones  with  bow  shock  interaction. 

3.  To  conduct  DNS  studies  of  nonlinear  wave  interaction  involving  Gortler  vortices 
in  hyi)ersonic  inlets  and  over  3-D  blunt  leading  edges  with  concave  surfaces. 

4.  To  conduct  the  DNS  studies  of  instability  and  transition  of  3-D  reacting  hyper- 
■sonic  boundary  layers  over  blunt  bodies.  The  fundamental  physics  of  the  real-gas 
('ffects  on  transition  are  investigated. 

5.  do  conduct  fundamental  studies  of  the  instability  properties  associated  with  hy¬ 
personic  boundary  layer  transition.  In  the  current  research  period,  the  focus  was 
mainly  on  the  linear  stability  analysis  of  hypersonic  flow  with  bow  shock  effects 
and  with  a  blunt  leading  edge. 
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3  ACHIEVEMENTS 


The  goal  of  this  research  project  is  to  develop  new  advanced  numerical  methods  and  to  per¬ 
form  DNS  studies  of  transient  hypersonic  reacting  flows  over  full  3-D  maneuvering  vehicles. 
The  DNS  tools  and  supporting  theoretical  approaches  are  used  to  gain  new  fundamental  un¬ 
derstanding  of  transition  phenomena  of  3-D  chemically-reacting  hypersonic  boundary  layers. 
Major  accomplishments  in  the  report  period  are  documented  in  the  publications  listed  in 
Section  7.  They  are  briefly  summarized  below. 

A.  DNS  of  Receptivity  of  Hypersonic  Flows  Over  Elliptical  Cones  and  Wedges. 

We  have  conducted  a  DNS  study  of  3-D  hypersonic  flows  over  elliptical  cross-section 
(  OIK'S  and  other  3-D  flows  by  using  our  14-node  IBM  SP-2  parallel  computer  purchased 
through  .‘XFOSR  funds.  We  have  conducted  DNS  studies  of  boundary  layer  receptivity  for 
•3-1)  hypersonic  flows  over  elliptical  cross-section  cones  and  over  blunt  leading  edges.  We 
are  currently  performing  parametric  DNS  studies  on  the  hypersonic  receptivity  physics  for 
li\  personic  flow  over  elliptical  blunt  cones  and  wedges.  Such  studies  are  valuable  in  prediction 
and  control  of  hypersonic  boundary  layer  transition. 

C.  DNS  of  Hypersonic  Gortler  Vortices  Over  Concave  Inlet  Surfaces 


We  liav('  finished  the  parametric  studies  of  the  linear  stability  analysis  of  Gortler  instabil¬ 
ity  in  h\  p('rsonic  boundary  layers.  Subsequently,  we  conducted  the  DNS  study  of  nonlinear 
(le\ clopiiK'nt  of  3-D  Gortler  vortices  on  concave  surfaces  with  the  effects  of  the  bow  shock. 
We  are  currently  pushing  the  simulation  domain  further  downstream  to  study  the  nonlinear 
breakdown  of  the  Gortler  vortices  in  hypersonic  boundary  layers. 

D.  DNS  of  Nonequilibrium  Reacting  Hypersonic  Flows 


We  studi('d  the  real-gas  effects  on  hypersonic  boundary-layer  flow  stability  and  transition 
l>y  using  the  D.NS  approach.  So  far.  we  have*  finished  the  reacting  flow  code  validation  and 
ha\'e  conduct('d  the  DNS  of  the  receptivity  of  2-D  reacting  hypersonic  boundary  layers  over 
l)lunt  l)o(lies.  We  are  currently  working  on  the  DNS  of  stability  and  nonlinear  breakdown 
of  iiy|)('rsonic  reacting  boundary  layers  over  a  flat  plate.  The  use  of  a  simpler  flat-plate 
geometry  enables  us  to  simulate  nonlinear  breakdown  better  and  to  gain  a  better  insight 
into  th(’  rc'al  gas  effects  on  the  reacting  hypersonic  boundary  layers. 

D.  Parallel  Implicit  Numerical  Algorithms 


During  the  report  period,  we  worked  on  developing  new  parallel  implicit  DNS  algorithms 
for  massively  parallel  computers  to  simulate  3-D  transient  reacting  hypersonic  boundary  layer 
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flows.  The  purpose  is  to  overcome  the  inter-node  communication  bottleneck  which  dramati¬ 
cally  reduces  parallel  simulation  efficiency  when  many  parallel  nodes  are  simultaneously  used 
in  simulations.  We  are  currently  investigating  the  possibility  of  using  a  divide-and-conquer 
parallel  algorithm  to  increase  the  parallel  efficiency  for  implicit  Navier-Stokes  solvers. 

E.  PSE  and  LST  Stability  Analysis  of  Non-Parallel  Hypersonic  Flows 

In  order  to  support  our  DNS  studies,  we  have  been  developing  PSE  and  LST  codes  for 
non-parallel  hypersonic  flow  with  bow  shock  effects.  So  far,  we  have  developed  and  validated 
two  independent  linear  stability  codes  and  have  extended  this  parallel  stability  tool  to  include 
the  non-parallel  effects  by  using  the  PSE  approach.  This  aspect  of  the  work  has  now  been 
finished. 

These  research  accomplishments  are  discussed  further  in  the  following  sections. 


4  ACCOMPLISHMENTS  AND  HIGHLIGHTS 


TIk'  major  accomplishments  and  highlights  of  our  research  are  summarized  in  this  section. 
Our  |)ubIications  in  the  reporting  period  are  listed  in  Section  7.  The  paper  numbers  cited  in 
this  s(‘ction  refer  to  those  used  in  Section  7. 

DNS  of  3-D  Hypersonic  Flow  Over  Elliptical  Cross-Section  Cones. 

Il>  ixusonic  flow  over  a  3-D  blunt  body  experiences  strong  lateral  pressure  gradients  that 
turn  the  streamlines  away  from  the  axial  direction,  inducing  a  skewed  boundary  layer  profile 
with  crossflow.  The  inflected  crossflow  velocity  profile  exhibits  an  inviscid  instability.  The 
(•hiiract('r  of  th(>  crossflow  instability  is  well-established  for  subsonic  and  moderately  super¬ 
sonic  Hows  over  swept  wings  and  certain  model  problems  such  as  rotating  disks,  cones,  and 
spheres  Stability  computations*'^'  have  shown  that  crossflow  increases  the  amplification 
rate  and  skc'wness  of  the  most  unstable  first  mode  wave,  and  the  most  amplified  second 
mode  wave  may  be  oblique  to  the  freestream  flow.  Recently,  Poggie  and  Kimmel  have 
doin'  ('X])crimental  work  to  examine  stability  and  transition  on  a  full  3-D  configuration  and 
to  demonstrate  a  case  in  which  the  crossflow  stability  had  a  significant  influence  on  boundary 
layer  transition. 

Our  purpose  is  to  conduct  numerical  studies  of  both  steady  and  transient  3-D  hypersonic 
flow  over  elliptical  cross-section  cones.  We  (Publications  [1,  11,  14])  have  extended  our 
previous  study  of  the  receptivity  of  hypersonic  boundary  layers  over  parabolic  blunt  leading 
edges  to  the  DNS  of  3-D  hypersonic  flows  over  the  cones.  The  first  test  case  is  the  receptivity 
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of  a  hypersonic  boundary  layer  to  freestream  monochromatic  planar  acoustic  disturbances 
for  a  Mach  15  and  Reynolds  number  6026  flow  over  a  2:1  elliptical  cross-section  blunt  cone 
as  shown  in  Fig.  1. 

Figure  2  shows  3-D  steady  solutions  for  Mach  number  contours  for  the  steady  viscous 
hypersonic  flow  over  a  2:1  elliptical  cross-section  blunt  cone  at  Mach  number  15.  The  bow 
shock  is  captured  exactly  as  the  outer  computational  boundary.  The  uneven  strength  of  the 
bow  shock  in  the  major  and  minor  axes  creates  the  circumferential  pressure  gradient  and 
cros.s  flow  velocity.  The  cross  mass  flow  thickens  the  boundary  layer  and  enlarges  the  standoff 
distance  in  the  minor  axis.  The  cross  flow  generates  a  streamwise  vorticity  component  in  the 
shock  layer.  Figure  3  shows  the  contours  of  streamwise  vorticity  components  in  a  x  ^constant 
grid  station.  The  figure  shows  that  the  vorticity  is  generated  by  cross  flow  velocity  on  the 
wall.  It  is  expected  that  the  cross  flow  instability  plays  an  important  role  in  the  transition 
(T  t  he  hypersonic  boundary  layers. 

Having  obtained  the  steady  flow  solutions,  the  generation  of  boundary-layer  w'aves  by 
freestream  acoustic  disturbances  is  simulated  for  3-D  hypersonic  flow  over  a  blunt  elliptical 
eoiK'  with  a  freestream  disturbance  wave.  The  receptivity  studies  are  carried  out  by  imposing 
acoustic  disturbances  in  the  freestream.  The  subsequent  interaction  with  the  bow  shock  and 
wa\(“  d('\('lopment  in  the  boundary  layer  are  simulated  by  using  the  full  nonlinear  Navier- 
■Stokes  ('(|ualions. 

I'igUK'  1  shows  the  3-D  contours  for  the  instantaneous  entropy  perturbation  s'  after 
the  uustead}’  computations  are  carried  out  for  enough  periods  in  time  that  periodic  solu¬ 
tions  ha\  ('  been  reached  in  the  entire  flow  field.  Figure  5  shows  instantaneous  entropy  and 
pressure  perturbations  in  the  major  and  minor  axes.  The  instantaneous  contours  show  the 
developiiH'ut  of  a  single  wave  mode  in  the  boundary  layer  along  the  surface.  The  receptivity 
simulation  r('sults  show  that  new  schemes  can  resolve  3-D  transient  hypersonic  flow  with 
physical  how-shock  oscillations  accurately. 

DNS  of  3-D  Hypersonic  Boundary  Layers  Over  Axisymmetric  Blunt  Cones 


In  addition,  we  have  also  extended  the  previous  planar  2-D  work  to  the  receptivity  of 
axisymmetric  boundary  layers  to  freestream  disturbances  in  Mach  15  flows  over  a  parabolic 
cone  (Publication  [17]).  The  receptivity  characteristics  of  axisymmetric  and  planar  hyper¬ 
sonic  boundary  layers  over  blunt  bodies  were  compared.  The  test  case  is  the  receptivity  of 
an  axis\  ininelric  boundary  layer  to  weak  freestream  acoustic  disturbance  waves  for  a  Mach 
15  hyp('rsonic  flow  past  a  parabolic  blunt  cone  at  zero  angle  of  attack.  The  axisymmetric 
steady  and  unsteady  solutions  for  velocity  vectors  are  shown  in  Fig.  6.  The  bow  shock  and 
develo|)ment  of  boundary  layers  along  the  body  surface  are  shown  in  this  figure.  The  figure 
also  shows  the  unsteady  instantaneous  perturbations  of  the  velocity  vectors  for  the  recep¬ 
tivity  of  the  axisymmetric  Mach  15  flow  over  the  parabolic  cone  with  zero  angle  of  attack. 
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The  results  show  that  new  schemes  can  resolve  3-D  transient  hypersonic  flow  with  physical 
bow-shock  oscillations  accurately.  The  receptivity  characteristics  of  axisymmetric  and  pla¬ 
nar  hypersonic  flow  over  blunt  bodies  are  studied  and  compared.  Compared  with  the  planar 
case,  the  axisymmetric  flow  over  a  blunt  cone  has  much  higher  after-shock  Mach  numbers 
and  much  fuller  boundary  layer  profiles.  Consequently,  the  axisymmetric  first  mode  wave 
generated  by  the  receptivity  process  has  longer  wave  lengths  and  smaller  growth  rates  than 
the  planar  case. 

DNS  of  Hypersonic  Gortler  Vortices  Over  Concave  Inlet  Surfaces 

The  concave  compression  surfaces  in  front  of  hypersonic  inlets  are  susceptible  to  Gortler 
(•(Hinter  rotating  instability  vortices,  in  addition  to  other  hypersonic  instability  waves  in 
hypersonic  boundary  layers  ^  \  Such  Gortler  instability  is  important  in  the  prediction  of  the 
iiiinsition  phenomena  in  inlet  flow.  In  addition,  the  Gortler  instability  can  also  be  exploited 
It)  ('iihance  the  instability  and  mixing  if  fuel  injection  is  applied  to  the  compression  surface 
bcfere  entrance  of  the  inlet.  The  DNS  studies  provide  an  ideal  tool  to  study  the  characteristic 
and  on-set  of  Gortler  vortices  for  hypersonic  flow  over  a  concave  surface  and  the  receptivity 
of  th<‘  Gbrtlcr  instability  to  the  disturbances  introduced  by  the  fuel  injection.  To  achieve 
inaxirinim  mixing,  it  is  possible  to  artificially  excite  the  Gortler  instability  by  specifying  the 
lucl  ini('ction  amplitudes  and  frequency.  We  have  been  working  on  the  DNS  of  hypersonic 
boundary  layer  stability  and  transition  over  a  concave  surface  with  the  effects  of  bow  shock 
intf'racti'on.  The  main  focus  is  to  gain  a  better  understanding  of  the  destabilizing  effects  of 
the  (!()rtler  vortices  on  the  boundary  layers.  In  this  project,  Gortler  vortices  in  hypersonic 
boundary  layers  are  investigated  mainly  by  DNS  while  LST  is  also  conducted  to  support  and 
help  validate  the  D.NS  analysis. 

In  the  K'port  period,  we  have  conducted  a  series  of  studies  on  linear  and  nonlinear 
interactions  of  Gortler  vortices  in  hypersonic  boundary  layers.  These  results  are  documented 
in  I’nblieaticnis  [>].  [0],  [12],  and  [13].  Some  sample  results  are  shown  in  this  section. 

We  have  obtained  both  steady  mean  flow  and  unsteady  flow  solutions  for  a  case  of  hy¬ 
personic  flow  over  a  blunt  concave  body,  where  the  mean  flow  is  two-dimensional  but  the 
utisteady  flow  of  the  Gortler  vortex  is  three-dimensional.  The  specific  flow  conditions  are 
.1/-,  =  15  and  Rcoc.  =  150753.175.  Both  linear  and  nonlinear  development  of  the  Gortler 
vort(‘X  in  the  boundary  layer  have  been  studied.  Figure  7  shows  steady  temperature  and 
pix'ssure  contours  for  Mach  15  mean  flow.  Both  figures  show  the  effects  of  concavity  of  the 
wall.  .After  obtaining  the  steady  flow  solutions,  DNS  of  the  linear  and  nonlinear  growth 
of  the  Gortler  vortices  are  simulated  by  imposing  inlet  disturbances  in  a  localized  compu¬ 
tational  domain  shown  in  Fig.  8.  The  inlet  disturbances  are  eigenfunctions  obtained  from 
linear  stability  analysis  using  the  simulated  mean  flow.  We  have  started  to  study  spatial 
Gortler  instability.  In  spatial  linear  stability  analysis,  eigenfunctions  do  not  depend  on  time, 
and  initial  disturbances  at  the  inlet  are  fixed  as  time  increases.  Initial  disturbances  prop- 


agate  spatially  and  converge  to  steady  state  condition.  In  the  current  computation,  the 
spanwise  nondimensional  wave  number  0  is  0.1,  which  is  nondimensionalized  by  boundary 
layer  thickness,  S.  Figure  8  shows  the  temperature  perturbation  contours  in  the  simulation. 
The  growth  of  the  Gortler  vortex  in  the  streamwise  direction  is  shown  by  the  increase  of  the 
intensity  of  the  disturbances.  The  numerical  solutions  agree  very  well  with  linear  predictions 
for  the  current  case  of  small  inlet  disturbances. 

In  addition,  we  have  also  developed  a  linear  stability  computer  code  to  study  the  Gortler 
instability  on  an  arbitrary  shape  surface.  VVe  have  validated  the  LST  code  for  the  Gortler 
instability  by  comparing  with  published  results  by  other  researchers  for  both  incompressible 
and  compressible  flows.  It  was  shown  that  the  code  gave  good  agreements  with  previous 
results.  Figure  9  shows  the  streamwise  velocity  contours  for  the  eigenfunctions  of  the  first 
and  second  Gortler  modes  for  an  incompressible  flow  case.  The  two  peaks  in  the  figure 
re|)resent  that  there  are  two  vorticity  layers  which  are  stacked  together.  Mode  II  contains 
one  more  layer  than  Mode  I.  Dashed  lines  represent  negative  values.  However,  higher  modes 
arc  rarrdy  observed  in  experiments  because  the  growth  process  of  these  modes  is  much  slower 
t  bail  one  of  the  primary  mode.  The  shapes  of  eigenfunctions  are  similar  to  those  in  Herbert 
I'i.uure  10  compares  eigenfunction  distributions  with  those  published  by  Kahawita  et  al  . 
flic  shape  of  eigenfunctions  obtained  by  our  calculations  agree  reasonably  well  with  those 
by  Kahawita. 

After  these  initial  studies,  we  have  extended  our  studies  to  the  nonlinear  development  of 
the  (iorth'f  \ortices  and  the  nonlinear  interaction  of  Gortler  vortices  with  other  instability 
modes,  file  nonlinear  breakdown  of  the  Gortler  vortrices  is  a  critical  step  in  the  laminar 
turbulent  transition  process  of  concave  hypersonic  boundary  layers.  Experiments  showed 
that  it  is  mainly  due  to  the  interaction  of  nonlinear  growth  of  Gortler  vortices  and  other 
forms  of  (list nrbatices.  In  the  nonlinear  growth  process,  mushroom  shaped  vortices  are 
liroduced  since  the  counter-rotating  vortices  pump  fluid  with  a  low  streamwise  velocity  away 
Irom  the  wall,  fhere  are  two  regions  in  Gortler  vortices  are  the  peak  region  (low  velocity 
region  I  atid  \alley  r('gion  (high  velocity  region).  These  two  regions  produce  the  mushroom 
shaped  \’ortices.  Interaction  between  these  vortices  and  traveling  waves  is  the  main  factor 
resulting  iti  biavik  down  to  turbulence. 

.As  an  example  of  the  DNS  results  on  nonlinear  development  of  Gortler  vortices  in  hy- 
per.sonic  bottndary  layers,  figure  11  and  12  show  streamwise  mean  velocity  distributions  as 
flow  moves  downstream.  The  development  of  mushroom  shaped  vortices  is  well  represented 
in  the  figure.  The  bow  shock  does  not  have  much  effect  on  the  flow  field,  since  Gortler 
vortices  develop  in  the  viscous  layers.  The  peak  and  valley  regions  are  clearly  shown  in  the 
(iguixu  While  the  middle  region  (peak)  tends  to  go  up,  others  become  narrower.  The  Gortler 
vortices  pump  vertically,  the  low-speed  fluid,  away  from  the  wall  in  the  peak  region  and  push 
the  high  speed  fluid  toward  the  wall  in  the  valley  region.  However,  there  is  the  limitation  of 
growitig  thickness  of  the  peak,  and  high  speed  fluid  starts  to  transfer  to  the  peak  region  near 
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the  W3.ll  3nd  the  low  speed  fluid  to  the  V3lley  3W3y  from  the  W3ll.  As  3  result,  the  mushroom 
shuped  vortices  3re  produced.  Detmls  of  the  results  C3n  be  found  in  [2],  [9],  [12],  3nd  [13]. 

DNS  of  Nonequilibrium  Reacting  Hypersonic  Flows 

For  high-enth3lpy  hypersonic  flows,  reul  g3s  effects  3re  expected  to  huve  3  strong  imp3ct 
on  the  flow  structure  3nd  he3ting  r3tes.  We  (Public3tions  [8]3nd  [15])  h3ve  been  working 
on  the  DNS  studies  of  inst3bility  3nd  tr3nsition  of  3-D  re3cting  hypersonic  boundury  kyers 
over  blunt  bodies  3nd  flut  pktes.  We  huve  extended  our  ide3l-g3S  DNS  shock-fitting  3-D 
codes  for  nonequilibrium  re3l-g3s  flow  simuktions.  The  effects  of  thermo-chemicul  re3Ctions 
on  bouiKkry-kyer  tr3nsition  3re  investiguted  by  both  DNS  3nd  theoreticul  3n3lysis  b3sed 
on  the  bound3ry  kyer  3ppro3ch.  The  re3l-g3s  model  is  3  five-species  3ir  model  with  both 
clieinieal  and  vibrational  nonequilibrium.  Figure  13  shows  steady  vorticity  contours  obtained 
b\  using  the  fifth-order  shock  fitting  khemes  for  Mach  15  hypersonic  reacting  flow  over  a 
blunt  parabolic  leading  edge.  The  real-gas  viscous  flow  solutions  are  compared  with  the 
corresponding  solutions  obtained  by  neglecting  the  real-gas  effects.  In  the  real  gas  case,  the 
bow  shock  is  more  significantly  curved  ,  and  the  standoff  distance  of  the  shock  is  smaller 
than  that  in  the  perfect  gas  case.  As  a  result,  larger  vorticity  is  generated  across  the  shock 
in  noii('(|uilibrium  flow.  The  figure  shows  that  the  real  gas  effects  have  a  strong  influence  on 
flow  propcMties.  The  real-gas  solution  has  much  stronger  entropy  and  vorticity  layers  which 
can  be  expected  to  induce  entropy  layer  instability  for  hypersonic  flows  . 

\\V  then  used  DNS  as  a  tool  in  studying  the  entropy  and  boundary  kyer  instability 
for  rca<  ting  3-D  hv'personic  flows  over  blunt  leading  edges.  In  the  unsteady  simulation, 
the  freest r<'am  disturbances  are  superimposed  on  the  steady  mean  flow  to  investigate  the 
de\  elopiiHuit  of  T-S  waves  in  the  boundary  kyer  with  the  effects  of  the  bow  shock  interaction. 
The  contours  of  vertical  velocity  perturbations  are  shown  in  Fig.  14. 

.Meanwhile,  we  have  also  finished  the  code  validation  and  evaluation  of  numerical  accuracy 
in  onr  hiiih-tuder  nonequilibrium  flow  code.  To  validate  this  numerical  method  and  computa- 
lional  code,  steady  hypersonic  flow  past  cylinders  at  the  same  flow  conditions  as  Hornung’s 
experiments  is  computed.  The  freestream  flow  conditions  are  partially  dissociated  nitrogen 
How  past  a  2  inch  diameter  cylinder  at  Ux  =  5590  in/s,  =  1833°  A,  Poo  =  2910  Pa,  92.7% 
.V.,  and  7.3/f  .V  by  mass,  T,^,  -  1833° A'.  Two  nonequilibrium  chemistry  models  were  used 
in  the  simulations  to  test  their  effects  of  real  gas  models  on  reacting  hypersonic  flows.  The 
first  moch'l  was  the  Park  model  and  the  second  model  was  the  Dunn-Kang  rate  coefficient 
model  Figure  15  compares  computational  results  using  the  two  models  (shown  in  the 
bottom  half  of  the  figure),  to  an  experimental  interferogram  published  by  Hornung  The 
comi)utational  shock  shape  and  standoff  distance  agree  very  well  with  the  experimental  data. 
The  inference  fringe  shapes  of  the  Park  Model  are  very  similar  to  the  experimental  photo.  For 
the  result  of  the  Dunn-Kang  Model,  there  are  visible  differences  between  the  inference  fringe 
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shapes  of  the  computation  and  experiment.  A  comparison  of  fringe  number  profile  along 
the  stagnation  line  between  the  two  chemical  models  is  shown  in  Figure  16.  Compared  with 
experimental  data,  the  Park  model  yielded  better  results  than  did  the  Dunn-Kang  model. 
The  results  indicated  that  the  choice  of  thermochemical  model  can  have  large  effects  on  the 
shock  shape  and  standoff  distance.  The  results  demonstrate  the  capability  of  our  numerical 
method  to  accurately  model  flows  in  thermochemical  nonequilibrium.  In  addition,  we  have 
done  extensive  validation  calculations  on  our  high-order  nonequilibrium  reacting  code,  in¬ 
cluding  the  comparison  of  stagnation  heating  rates  for  both  perfect-gas  and  nonequilibrium 
flows  with  the  classical  theory  of  Fay  and  Riddell  The  theoretical  prediction  in  the  heat 
transfer  at  the  stagnation  point  agrees  well  with  numerical  simulations  of  full  Navier-Stokes 
('(|uations. 

riu'  real  gas  effects  on  hypersonic  boundary  layer  stability  and  transition  involve  com- 
pl('.\  flow  physics  of  nonequilibrium  reaction  and  relaxations.  In  order  to  better  study  the 
inecliaiiism  of  real  gas  effects  on  boundary  layer  stability  and  transition,  we  have  focused 
oil  a  simple  flow  domain  of  a  flat  plate  with  nonequilibrium.  Again,  DNS  is  used  to  study 
tin-  liiK'ar  and  nonliner  stabilities  of  nonequilibrium  reacting  flow  over  a  flat  plate.  Some 
initial  |■(‘sults  have  been  presented  in  Publication  [15].  The  numerical  solutions  are  validated 
by  l)()nndarv  layer  theory  and  linear  stability  analysis.  Figure  17  shows  the  contours  of 
instantaneous  pressure  perturbation  after  the  flow  field  reaches  a  periodic  state  for  a  case  of 
dexclopment  of  instability  waves  in  a  flat  plate  boundary  layer. 

Parallel  Implicit  Numerical  Algorithms  for  DNS  of  Reacting  Hypersonic  Flows 


In  Publications  [5],  [10],  [20]  and  [21],  we  continue  to  develop  and  validate  new  parallel 
computer  code's  using  our  new  semi-implicit  schemes  for  the  DNS  of  3-D  hypersonic  boundary 
laye'is  ()V('r  l)lunt  bodies. 

In  [lO;.  we  (h'veloped  new  parallel  implicit  D.NS  algorithms  for  massively  parallel  com¬ 
puter'-  to  simulate  ;5-l)  transient  reacting  hypersonic  boundary  flows.  The  computational 
domain  is  (lecom|)osed  into  a  number  of  smaller  computational  domains.  Parallel  program¬ 
ming  with  the  .Message-Passing  Interface  (.MPI)  model  is  applied  to  parallelize  the  3-D  code 
to  ((unpute  the  different  domains  by  using  different  computer  processors  respectively  at  the 
sanu'  time.  We  used  a  divide-and-conquer  parallel  algorithm  to  increase  the  parallel  efficiency 
for  implicit  .Xavier-Stokes  solvers.  The  purpose  is  to  overcome  the  inter-node  communica¬ 
tion  bottleneck  which  dramatically  reduces  parallel  simulation  efficiency  when  many  parallel 
.nodes  ar<'  simultaneously  used  in  simulations.  From  the  results  of  studying  the  convection- 
diffusion  model  equation  and  the  instability  of  supersonic  Couette  flow,  the  computational 
('fliciency  can  be  improved  by  using  parallelized  semi-implicit  methods  while  maintaining 
the  same  accuracy  as  that  of  the  parallelized  explicit  methods.  We  can  also  get  the  satisfac¬ 
tory  results  by  using  the  new  parallelized  semi-implicit  methods  for  the  direct  simulations 
of  hypersonic  boundary  layers  over  blunt  bodies  to  freestream  acoustic  disturbances. 
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During  the  report  period,  a  set  of  new  high-order  upwind  schemes  with  and  without  the 
shock-fitting  procedure  were  also  validated  and  tested  in  computations  for  the  linear  wave 
equation  and  nonlinear  Navier-Stokes  equations.  The  numerical  tests  show  that  the  new 
schemes  can  achieve  efficient  and  high-order  numerical  simulation  for  reacting  hypersonic 
boundary  layers  over  three-dimensional  blunt  bodies.  A  semi-implicit  high-order  upwind  fi¬ 
nite  difference  shock  fitting  method  has  been  tested  and  validated  in  Ref.  [12],  which  includes 
several  test  cases:  TD  linearized  model  convection-diffusion  equation,  2-D  DNS  of  stability 
of  supersonic  Couette  flow,  2-D  DNS  of  stability  of  supersonic  flat  plate  boundary  layer, 
and  the  DNS  of  receptivity  to  freestream  acoustic  disturbances  for  2-D  hypersonic  boundary 
layer  over  a  parabola.  The  CPU  time  reduction  for  the  semi-implicit  method  for  both  cases 
is  typically  in  the  range  of  8-10  times  over  those  required  by  the  explicit  calculations  with 
the  same  numerical  accuracy. 

\\'(>  have  also  conducted  code  validation  on  hypersonic  flow  over  an  axisymmetric  body 
at  zero  angle  of  attack  for  the  new  high-order  parallelized  shock-fitting  code.  A  sphere  is 
chosi'ii  as  the  validation  case  since  there  are  a  lot  of  experimental  results  for  a  sphere  in 
supi'isonic  flow  Grid  refinement  comparisons  are  also  used  to  ensure  that  solutions 

are  grid  independent.  The  flow  conditions  are  Moo  =  5.25  and  Rej,  =  36,159.3.  Figure  20 
shows  the  steady  solution  for  a  set  of  90  x  60  x  32  computational  grids  obtained  by  the 
numerical  simulations  for  Mach  number  5.25  flow  over  a  sphere.  Since  these  steady  flows  are 
axisymmetric.  we  present  the  results  only  in  a  computational  surface  of  a  fixed  azimuthal 
angif'.  I  his  figure  also  shows  the  contours  of  steady  axisymmetric  flow  solutions  of  Mach 
numbers  and  the  computed  and  experimental  pressure  coefficients.  The  agreement  is  good 
within  th<‘  scatter  in  the  experimental  data. 

PSE  Studies  of  Non-Parallel  Hypersonic  Flow  Instability 

III  order  to  support  our  major  efforts  in  the  DNS  studies,  we  have  developed  two  linear 
siabilitv  code's  for  hypersonic  boundary  layer  shear  flows  and  have  been  conducting  LST 
studies  of  h\])er.sonic  wave  characteristics.  In  Publications  [16]  and  [19],  we  studied  the 
lineiir  stabilitv  of  3-D  hypersonic  boundary  layers  over  blunt  leading  edges.  The  linear 
stability  ainilysis  is  performed  using  a  global  spectral  collocation  method  accounting  for  the 
shock  effects  by  using  Rankine-Hugoniot  shock  conditions  on  the  upper  boundary.  It  is 
shown  that  iti  addition  to  the  boundary  layer  first  modes  and  higher  modes  (Mack  modes), 
the  liiu-ar  stability  of  the  hypersonic  flow  between  a  bow  shock  and  a  parabolic  leading  edge 
has  a  iK'w  family  of  modes,  shock  modes.  The  shock  modes  are  important  mainly  in  the 
shock  la\er.  In  addition,  the  boundary  layer  modes  are  also  affected  by  the  existence  of 
the  shock.  In  comparison  with  the  DNS  results,  it  is  found  that  both  the  LST  and  DNS 
resoh'e  the  boundary  layer  modes  well  in  terms  of  the  wave  number  and  the  shape  function 
of  the  pressure  disturbance.  However,  disagreement  in  the  growth  rates  between  LST  and 
DNS  solutions  has  been  observed.  The  difference  may  be  due  to  the  fact  that  the  underlining 
assumptions  of  LST  is  the  parallel  assumption,  which  may  not  be  proper  for  the  blunt  bodies. 
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PSE  (parabolized  stability  equations)  methods,  which  account  for  weak  non-parallel  eflfects, 
are  currently  being  developed  to  improve  the  comparisons  with  DNS  for  hypersonic  blunt 
body  flows.  The  purpose  of  using  PSE  in  addition  to  DNS  is  to  provide  an  alternative  means 
of  analyzing  the  stability  of  the  hypersonic  blunt  body  flows  where  nonparallel  effects  are 
important. 

We  have  since  developed  a  nonparallel  PSE  code  for  hypersonic  flow  over  blunt  bodies. 
C'urrently,  some  benchmark  tests  have  been  conducted  for  a  flat  plate  supersonic  boundary 
layer.  The  results  shown  in  Fig. 18  are  the  growth  rates  vs  R  at  F  =  40,h  =  0.15,  To  =  311/v' 
where  b  =  /.?*  10^//?o,  Ro  is  the  initial  Reynolds  number.  Our  PSE  results  differ  slightly  from 
the  results  from  Bertolotti  which  may  be  due  the  use  of  different  viscosity  models  and 
different  approaches  in  obtaining  initial  conditions.  Nonetheless,  our  PSE  results  agree  with 
th('  result  from  El-Hady’s  multiple  scale  analysis  Another  case  with  6  =  0  is  computed 
with  both  LST  and  nonparallel  PSE,  wavenumbers  and  growth  rates  from  both  methods 
are  shown  in  left  and  right  sides  of  Fig.  19  respectively.  Results  from  Bertolotti  are  also 
sIkavii  to  compare  with  our  growth  rate  data.  Notice  that  Bertolotti’s  results  with  parallel 
l)asic  flow  are  LST  results.  The  difference  again  is  due  to  the  reasons  mentioned  earlier  as 
one  can  s('e  that  the  use  of  Pr  =  0.72  gives  a  better  match  with  Bertolotti’s  data  where  Pr 
is  not  a  constant,  than  that  of  Pr  =  0.7.  The  nonparallel  effects  in  this  case  stabilize  the 
Ihnv  as  the  growth  rates  from  PSE  are  observed  to  be  smaller  than  the  LST  results.  The 
wavcMiumbers  from  both  methods  are  quite  close. 


5  FUTURE  WORK 


So  far.  w('  liave  developed  new  methodologies  and  have  demonstrated  the  feasibility  of  our 
approach  as  a  |)ractical  DNS  tool  to  explore  transient  hypersonic  reacting  flow  physics. 
Th('  focus  of  future  work  is  on  the  DNS  studies  the  hypersonic  flow  physics  related  to  the 
h\  |)ersouic  Itonndary  layer  transition,  and  to  extend  these  DNS  and  theoretical  studies  to  3- 
I)  transic'iit  hy|)ersonic  flows  with  complex  real-gas  models  and  complex  geometries,  and  the 
D.\S  of  nonlinear  breakdown  and  interactions  of  instability  modes  in  hypersonic  boundary 
layers.  The  spc'cific  tasks  are: 

1.  'fo  perform  DNS  of  the  receptivity  and  transition  of  3-D  hypersonic  boundary 
layers  over  blunt  elliptical  cross-section  cones  with  cross  flow  instability,  nonlinear 
l)reakdown.  and  general  3-D  geometries. 

2.  'I'o  perform  DNS  of  the  receptivity  and  transition  of  3-D  nonequilibrium  reacting 
hypersonic  boundary  layers  over  flat  plate,  blunt  leading  edges,  and  blunt  cones, 
d'he  real  gas  effects  on  linear  development  and  nonlinear  breakdown  of  instability 
modes  will  be  studied  by  the  DNS  approach. 
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3.  To  perform  DNS  of  the  receptivity  and  transition  of  hypersonic  boundary  layers 
over  concave  blunt  bodies  with  Gortler  instability,  as  a  first  step  in  the  DNS  of 
complete  inlet  flow  fields  with  fuel  injection.  The  focus  is  on  the  late  stage  of 
nonlinear  breakdown  and  interaction  of  the  Gortler  and  other  instability  modes. 

4.  To  develop  new  accurate  and  efficient  parallel  numerical  algorithms  for  reacting 
hypersonic  flow  DNS  schemes.  The  main  focus  is  on  the  development  of  efflcient 
parallel  implicit  algorithms  for  massively  parallelized  computers. 
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10  FIGURES 


Figure  1:  Hypersonic  flow  flelcl  in  the  direct  numerical  simulation  of  3-D  reacting  hypersonic 
boundary-layer  transition  and  receptivity  to  free-stream  disturbances  over  a  blunt  elliptic 
cross-section  cone. 
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I*it»urc  1:  Instantaneous  3-D  entropy  perturbations  induced  by  freestream  planar  acoustic 


di'>t  urbanc('s. 
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Figur('  o:  Instantaneous  pressure  perturbations  induced  by  freestream  planar  acoustic  dis¬ 
turbance's  in  the  major  and  minor  axes. 


Figurt'  (i:  Solutions  of  steady  and  unsteady  simulations  of  axisymmetric  Mach  15  viscous 
hypersonic  flow  over  a  blunt  cone  (Upper  figure:  steady  velocity  vectors,  lower  figure:  in¬ 
stantaneous  perturbations  of  velocity  vectors  caused  by  the  receptivity  to  freestream  sound 
at  a  nondimensional  frequency  F  =  1770). 
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l  iguK*  7:  l<'mj)rratiire  and  pressure  contours  for  the  steady  hypersonic  mean  flow  over  a 
concav.-  blunt  body  at  =  lo.  7’;  =  101.059  K,  P;^  =  10.3  Pa,  T’  =  1000  K.  Re^  = 
15075.3.17. 
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Coinpul  at  ional 
Domain 


iourr  8-  A  schematic  of  the  concaved  body(engine  inlet)  on  which  shock  and  Gortler 
art  ices  exist.  Temperature  contours  from  DNS  results,  showing  Gortler  vortex  growth 
cer  a  concaved  bodv  surface  on  which  shock  and  Gbrtler  vortices  exist.  Flow  conditions  are 


=  15.  T. 


101.059  K.  p: 


10.3  Pa,  r;  =  1000  K,  Re^  =  150753.17 
23 


I'igure  10;  Eigenfunction  distributions  of  the  Gortler  instability  first  mode,  for  incompressible 
flow  over  a  concave  wall,  in  two  different  sets  of  grids  at  G  =  2.569,  R  =  3.2m,  =  292.4/tr, 

r,;.  =  'ymiscc,  and  j3  —  1.0  (  Streamwise  velocity,  normal  velocity,  spanwise  velocity,  and 
pressure).  Results  are  compared  with  Kahawita  et  al(1979). 


I'iuurf'  1  1:  Distributions  of  iso-contours  of  streamwise  mean  velocity  along  the  streamwise 
(lir(H  tit)n  for  /one  7.  8  and  9.  The  size  of  grids  is  483  x  121  x  64.  Mushroom  shaped  vortices 
dr\(‘|(>i).  as  flow  moves  downstream,  due  to  nonlinear  effects. 


Figure'  12:  Distributions  of  iso-contours  of  streamwise  mean  velocity  along  the  streamwise 
dire'ction  for  zone  10  and  11. 
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Top  Half:  Perfect  Gas. 
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Bottom  Half:  Real  Gas. 


I'igurc  l;l:  \orticity  contours  of  2-D  hypersonic  (Moo  =  15)  viscous  flow  over  a  parabolic 
leading  ('(lg('.  where  the  reacting  flow  solution  is  compared  with  the  perfect  gas  solution.  The 
solutions  are  obtained  by  using  fifth-order  shock  fitting  codes  with  nonequilibrium  source 
terms.  The  real  gas  effects  introduce  much  stronger  vorticity  in  the  entropy  shear  layer 
behind  the  curve  bow  shock. 
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l'igur('  I  I:  rnsteady  vertical  \Tlority  contours  of  2-D  hypersonic  {Moo  —  ho)  viscous  flow 
ov(M’  a  parabolic  leading  edge.  The  solutions  are  obtained  by  using  fifth-order  shock  fitting 
cod('s  with  nonecjuilibriuni  source  terms. 


1*  idlin'  15:  C’omparison  of  computed  pressure  contours  (bottom)  with  experimental  results 
of  lloriiung  (1991)  (top),  (a)  Park  model,  (b)  Dunn-Kang  model  for  hypersonic  flow  over  a 
cylind('r. 
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Figure  16:  Computed  fringe  numbers  for  the  two  chemistry  models  compared  with  experi¬ 
mental  values  of  Hornung  (1972)  along  the  stagnation  line  for  hypersonic  flow  over  a  cylinder. 
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17:  Contours  of  instantaneous  pressure  perturbation  for  the  case  of  uj/loq  = 


HB,  PSE 
El-Hady 
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Fi-iirr  18:  Growth  rates  vs  R  at  F  =  40,6  =  0.15,  To  =  311/\  at  Mach  1.5  computed  by 
I’Sl-:  for  flow  over  a  flat  plate.  The  reference  length  scale  is  the  local  scale  5.  Nonparallel 
l)a.sic  How  results. 
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Figure  19:  Growth  rates  (right)  and  wavenumber  (left)  from  LST  and  nonparallel  PSE  vs  R 
at  F  =  40,6  =  0,ro  =  311/\  at  Mach  1.5,  flat  plate  boundary  layer.  The  reference  length 
scale  is  the  local  scale  5.  Results  from  Bertolotti  are  also  shown  in  symbols. 


Figun'  20:  St('a(ly  flow  solutions  for  M~^  =  5.25  axisymmetric  steady  flow  over  a  sphere 
In'  using  high  order  3-1)  shock-fitting  method,  (Upper  figure:  grids,  middle  figure:  Mach 
numl)(u'  contours,  and  lowc'r  figure:  pressure  coefficients  on  body  surface). 
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